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Abstract 

In this paper, we consider random access by nodes that have energy harvesting capability. Each node is equipped 
with both a queue for storing the arriving packets and a battery for storing the harvested energy chunks, where the 
packet arrival and the energy harvesting events are all modeled as discrete-time stochastic processes. In each time 
slot, each node attempts to transmit the head-of-the-line packet in the queue with some probability if its battery 
is non-empty, and each transmission consumes one chunk of energy. Therefore, the transmission by one node is 
not just limited by the availability of packets in the queue, but by the availability of energy chunks in the battery. 
In most of related previous work, it was implicitly assumed that there exists unlimited energy for transmission, 
which is impractical in many distributed systems such as the ones based on sensors. In this work, we characterize 
the exact stability region when a pair of bursty nodes, which are harvesting energy from the environment, are 
randomly accessing a common receiver. The analysis takes into account the compound effects of multi-packet 
reception capability at the receiver. The contributions in the paper are twofold: first, we accurately assess the effect 
of limited, but renewable, energy availability due to harvesting on the stability region by comparing against the case 
of having unlimited energy. Secondly, the impact of the finite capacity batteries on the achieved stability region is 
also quantified. 

Index Terms 

Stochastic energy harvesting, stability, interacting queues, multipacket reception capacity, random access 



I. Introduction 

EXPLOITING renewable energy resources from the environment, often termed energy harvesting, permits 
unattended operation of infrastructureless distributed wireless networks. There are various forms of energy 
that can be harvested including thermal, vibration, solar, acoustic, wind, and even ambient radio power 0]]- 
0. These aspects of energy harvesting are not examined here. Despite the rapid advancement of hardware 
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technologies, including devices such as solar cells, thermoelectric generators, and piezoelectric actuators, the study 
of communication systems comprised of nodes that have energy harvesting capability are still in a very early 
stage. In 01, the capacity of the additive white Gaussian noise channel with stochastic energy harvesting at the 
source node was shown to be equal to the capacity of the channel with an average power constraint equal to 
the energy harvesting rate. Like most information-theoretic research, however, the result is obtained for a single 
source-to-destination channel with an infinitely backlogged source and is asymptotic in nature. 

When dealing with nodes equipped with non-rechargeable batteries, the common objectives were usually short- 
term such as maximizing the finite network lifetime 121, 0. The additional functionality of harvesting energy permits 
our assessment of the system long-term performance such as throughput, fairness and stability. On the other hand, 
the ALOHA protocol, the simple scheme of attempting transmission randomly, independently, distributively, and 
based on simple ACK/NACK feedback from the receiver, has gained continued popularity since its creation by 
Abramson |7]. It is especially suitable for distributed multi-access communication systems due to its simplicity 
and the independence of the centralized controller. It also serves as a cornerstone benchmark for assessment of 
performance of more elaborate schemes. In such a context, we revisit the canonical problem of the random access 
stability when nodes are powered by batteries recharging from randomly time-varying renewable energy sources. 
To this end, we characterize the stability region of the system which is defined as the set of packet arrival rate 
vectors for which all queues in the system are stable 3- 

In order to put our contribution in perspective, we start with some background on the stability of random access 
systems. The characterization of the stability region of random access systems for bursty sources (in contrast to 
infinitely backlogged sources, for which the concept does not make sense) is known to be extremely difficult even 
with unlimited energy for transmissions. This is because each node transmits and thereby interferes with the others 
only when its queue is non-empty. Such queues are said to be interacting, or coupled, with each other in the sense 
that the service process of one depends on the status of the others. Consequently, the individual departure rates of the 
queues cannot be computed separately without knowing the stationary probability of the joint queue length process, 
which is intractable [8]. This is the reason why most work has focused on small-sized networks and only bounds or 
approximations are known for the networks with larger number of nodes ll8l- lfT3l . In |9), the exact stability region 
was obtained for the two-node case and for an arbitrary number of nodes with symmetric parameters (that is, equal 
arrival rates and random access probabilities). In |8 ], a sufficient stability condition for an arbitrary number of nodes 
with asymmetric parameters was obtained. In [10], the necessary and sufficient stability condition was derived but it 
can only be evaluated up to the three-node case. The concept of the instability rank was introduced in ifTTl to further 
improve the inner bound for the general asymmetric cases. In lfT2l . an approximate stability region was obtained 
for an arbitrary number of nodes based on the mean-field asymptotics. All the above results were derived under the 
collision channel model in which, if more than one nodes transmit simultaneously, none of them are successful. This 

'A queue is said to be stable if it reaches a steady state and does not drift to infinity. A formal definition is given in Section ITT1 
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is too pessimistic assumption today in the sense that a transmission may succeed even in the presence of interference 
|[T3l - |[T5l . In such a context, the two-node stability result was extended to the channel with the multipacket reception 
(MPR) capability which enables the probabilistic reception of simultaneously transmitted packets [13]. On the other 
hand, all existing work on the stability of random access systems implicitly assumes unlimited availability of energy 
for transmissions, which is impractical for several distributed communication systems. 

In this work, we focus attention on the effect of limited availability of energy in each node's battery, that can 
be recharged by harvesting energy from the environment, on the system stability region when a pair of such nodes 
are randomly accessing a common receiver with MPR capability. Note that the analysis becomes significantly more 
challenging than in previous approaches, because the service process of a node depends not only on the status of 
its own queue and its battery, but also on the status of the other node's queue and battery. The key fact that makes 
the analysis tractable in this doubly interacting system is that the energy consumption is somewhat simplified as it 
does not depend on the success of the corresponding transmission. For the characterization of the stability region, 
we first obtain an inner and an outer bound of the stability region for a given transmission probability vector. 
Since an input rate vector that is outside of the stability region at a given transmission probability vector may be 
stably supported by another transmission probability vector, determination of the closure of the stability region is 
necessary and important. Consequently, we take the closure of the inner and the outer bound separately over all 
feasible transmission probability vectors. The remarkable result is that they turn out to be identical. Therefore, our 
characterization is exact in the sense that the bounds are tight in terms of the closure. It needs to be pointed out 
that even without the energy availability constraints, the exact characterization of the stability region in terms of the 
closure is known only for the two-node case, as also considered here. Finally, we remark that the results presented 
in this work generalize those of previous work that assumed unlimited energy for transmission over the collision 
channel flS), and over the channel with MPR capability |fT3l . 

The rest of the paper is organized as follows. In Section [TTJ we revisit the notion of the stability, describe the 
communication protocol, and explain the channel model and the packet arrival and energy harvesting processes. In 
Section [Till we present our main result on the stability of the slotted ALOHA with stochastic energy harvesting 
when nodes are equipped with infinite capacity batteries. The proof of our main result is delayed until Section [IVJ 
In Section |VJ we consider the case when nodes are equipped with finite capacity batteries and identify the impact 
of that restriction on the stability region. Finally, we draw some conclusions in Section [VT] 

II. System Model 

A. General Description of the Model 

The system consists of a pair of source nodes randomly accessing a common receiver, each powered by its own 
battery that can recharge from randomly time-varying renewable energy sources as shown in Fig. Q] Each node 
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Fig. 1. Random access of nodes with stochastic energy harvesting 
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has an infinite size queue for storing the arriving packets, that have fixed length, and a batter)o for storing the 
harvested energy. Time is slotted and the slot duration is equal to one packet transmission time. Energy is harvested 
in chunks of fixed size and one chunk of energy is consumed in each transmission. That is, the size of the chunk 
is equal to the slot duration times the power needed to transmit the fixed size packet over the slot duration. Let 
(-Aj(n), n > 0) and (Hi(n),n > 0) denote the packet arrival and energy harvesting processes at node i, respectively. 
They are modeled as independent and identically distributed (i.i.d.) Bernoulli processes with E[Ai(n)] = Aj and 
E[Hi(n)] = 5i. The processes at different nodes are also assumed to be independent of each other. Let Qi(n) and 
Bi (n) represent the number of buffered packets and the number of stored energy chunks at node i at the beginning 
of the n-th slot, respectively. Then, Qi(n) and Bi{n) evolve according to 



and 



Qi(n + 1) = Qi{n) - m(n) + Ai(n) 



Bi(n + 1) = Bi(n) - L(n) + H t {n) 



(1) 



(2) 



where /ij(n) £ {0, 1 } is the actual number of packets that are successfully serviced and, thus, depart from the 
queue of node i during time slot n, and L(n) is the indicator function such that L(ra) = 1 if node i transmits 
at time slot n, while otherwise L(ra) = 0. Node i is said to be active if both its data queue and its battery are 
non-empty at the same time, so that it can then transmit with probability p^. If either the queue or the battery is 
empty, node i is said to be idle and remains silent. Since the nodes are accessing a common receiver, the service 
variable /J,i(n) depends not only on the status of its own queue and of its battery but also on the status of the other 



2 The battery capacity is first assumed to be infinite and later relaxed to any finite number. 
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node's queue and battery. Furthermore, it also depends on the underlying channel model. The collision channel 
model, for instance, excludes all the possibility of success when packets are transmitted simultaneously. 

B. Multipacket Reception Channel 

The channel model with MPR capability considered in this work is a more realistic and general form of a packet 
erasure model which captures the effect of fading, attenuation and interference at the physical layer, along with 
the capability of multi-user detectors at the receiver |fT3l - lfT6l . Denote with g^i^ the success probability of node 
% when a set M. of nodes are transmitting simultaneously. The success probability qi\M is related to the physical 
layer parameters through 

Qi\M = Pr [7i|A4 > 0] (3) 

where 7ji^ denotes node i's signal-to-interference -plus-noise-ratio (SINR) at the receiver given set Ai of transmit- 
ters and 9 is the threshold for the successful decoding of the received packets, which depends on the modulation 
scheme, target bit-error-rate, and the number of bits in the packet, i.e., the transmission rate. Once the channel 
statistics are known, the packet reception probabilities can be readily computed. In Appendix [A] the MPR proba- 
bilities are obtained in a Rayleigh fading environment. Of course, equation ([3]) is an approximation since it does 
treat interference as Gaussian noise; however, it is used widely and represents a compromise between accuracy and 
cross-layer modeling [17]. 

C. Stability Criteria 

We adopt the notion of stability used in ifTOl in which the stability of a queue is equivalent to the existence of 
a proper limiting distribution. In other words, a queue is said to be stable if 

lim Px[Qi(n) < x] = F(x) and lim F(x) = 1 (4) 

n— )-oo x— >oo 

If a weaker condition holds, namely, 

lim liminf Pr[Qj(n) < x] = 1 (5) 

x— >oo n— >oo 

the queue is said to be substable or bounded in probability. Otherwise, the queue is unstable. A stable queue is 
necessarily substable, but a substable queue is stable if the distribution tends to a limit. If Qi(n) is an aperiodic 
and irreducible Markov chain defined on a countable space, which is the case considered in this paper, substability 
is equivalent to the stability and it can be understood as the recurrence of the chain. Both the positive and null 
recurrence imply stability because a limiting distribution exists for both cases although the latter may be degenerate. 
Loynes' theorem, as it relates to stability, plays a central role in our approach [18]. It states that if the arrival and 
service processes of a queue are strictly jointly stationary and the average arrival rate is less than the average 
service rate, the queue is stable. If the average arrival rate is greater than the average service rate, the queue is 
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unstable and the value of Qi(n) approaches infinity almost surely. If they are equal, the queue can be either stable 
or substable but in our case the distinction is irrelevant as mentioned earlier. 

At a given energy harvesting rate and for given transmission probability vectors, the stability region 6(5, p) is 
defined as the set of arrival rate vectors A = (Ai, A2) for which all queues in the system are stable. The stability 
region of the system &(S) is defined as the closure of 6(6, p) over all possible transmission probability vectors, 
i.e., 6(5)±{J pem2 6(6,p). 

III. Main Result 

This section presents our main results on the stability of slotted ALOHA for the two-node case with stochastic 
energy harvesting when the capacity of the batteries is assumed to be infinite. Define Aj = qi\{i} — ?j|{i,2}> which 
is the difference between the success probabilities when node i is transmitting alone and when it transmits along 
with the other node j (7^ i). The quantity Aj is strictly positive since interference only reduces the probability of 
success. Let us define the following points in the two-dimensional Euclidean space to facilitate the description of 
our main theorem: 

PA = (0,5 2 q 2 \{2}) (6) 
P B = ( g2|{2}(gi|{i} ~ Ai<5 2 ) 2 Ai£|g 2 | {2} ^ 




Pb 2 =\ — —) (8) 

Q2\{2} ^1<?2|{2} 

Pb 3 = {h(qi\{i} ~ Aifc), h(q 2 \{2} ~ A 2 5i)) (9) 

Pc= (5i?i|{i},0) (10) 

where Pb 1 , Pb 2 , and Pb 3 are in the first quadrant and Pa and Pq are on y and x-axes, respectively. These points 
can be seen in Fig. [2] Let us further define 

+ M (ID 

9l|{l} Q2\{2} 

which is non-negative and decreasing as the MPR capability improves. 

Theorem 3.1: If ^ > 1, the boundary of the stability region &(S) of the slotted ALOHA at a given energy 
harvesting rate S is described by three segments: (i) the straight line connecting Pa and Pb x , (ii) the curve 



VA 2 Ai + V A 1 A 2 = ^/<?1|{1}<?2|{2} (12) 

from Pb 1 to Pb 2 , and (iii) the straight line connecting Pb 2 and Pq- If ^ < 1, it is described by two straight lines: 
(i) the line connecting Pa and Pb 3 and (ii) the line connecting Pb 3 and Pq- 

Proof: The proof is presented in the next section. ■ 
In Fig. 12 we illustrate the stability region &(S) for different packet reception probabilities. The boundary of the 
region is indicated by the solid line. The case with unlimited energy, i.e., 5i = 1, Vi G { 1, 2 }, is also depicted in 
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Fig. 2. Two-node stability region &(S) at different MPR probabilities where <5 = (0.8,0.7) and qi\{iy = 0.9, q 2 |{i} = 0.8 
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the figure with the dotted line. The difference between the two regions, therefore, can be understood as the loss due 
to the limited availability of energy imposed by the variable battery content and the stochastic recharging process. 
Corollary 3.1: If ^ > 1, the stability region &(S) is non-convex, whereas if * < 1, it is a convex polygon. 
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When ^ = 1, the region becomes a right triangle. 

This corollary can be easily verified by comparing the slopes of the lines from Pa to Pb 1 and from Pb 2 to Pc 
and those from Pa to Pg 3 and from Pg 3 to Pc*. Specifically, when \& = 1, the curve (fT2l ) shrinks to a point whose 
coordinates are identical to those of Pb 1 and Pb 2 and the slopes of the lines from Pa to Pb x and from Pb 2 to Pc* 
become identical. 

When we do not have MPR, i.e., when the channel is described by the classical collision channel model, we can 
obtain the stable region as described in the following corollary and shown in Fig. [3] 

Corollary 3.2: The stability region &(5) of the slotted ALOHA under the collision channel model is described 
as follows. If 5± + 82 > 1, its boundary is described by three segments: (i) the line segment connecting (0, £2) and 
((1 - 5 2 ) 2 , <5|), (ii) the curve ^/\[ + y/X^ = 1 from ((1 - 5 2 ) 2 ,S 2 ) to (<5f,(l - 5i) 2 ), and (iii) the line segment 
connecting (5 2 , (1 — 5i) 2 ) and (<5i,Q). If Si + 62 < 1, it is described by two lines: (i) the line segment connecting 
(0,6 2 ) and (<5i(l — 82), 62(1 — Si)) and (ii) the line segment connecting (<5i (1 — 62), 62(1 — 61)) and (<5i,0). 
The corollary is obtained by substituting ftim = 1 and qi\{i t 2} = 0, Vi € { 1, 2 }, into Theorem 13. II 

Corollary 3.3: The stability region 6(1) of the slotted ALOHA with the unlimited energy for transmission under 
the collision channel is the region below the curve \/Ai + \/^2 = 1 in the first quadrant of the two-dimensional 
Euclidean space. 

This last corollary, which is obtained by substituting Si = 1, Vi € {1,2}, into Corollary 13.21 reconfirms the 
well-known result on the stability of the slotted ALOHA obtained in (H. 
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Fig. 4. An inner and an outer bound for <S(S,p) 



IV. Stability Analysis 

In this section, we prove our main result presented in the previous section. We first derive a sufficient condition 
for stability in Section IIV-AI and, separately, a necessary condition for stability in Section IIV-BI for given energy 
harvesting rates S and transmission probabilities p, which yield an inner and an outer bound of &(6, p), respectively, 
and they are shown in Fig. |4] The achievability and the converse of Theorem 13.11 is shown in Section IIV-CI by 
taking the closure of the inner and the outer bounds of &(S,p) over p and by observing that these closures turn 
out to be identical. 

A. Sufficient Condition 

For the sufficiency, we show that any arrival rate vector A = (Ai , A2) that is componentwise less than the saturated 
throughput vector of the system, denoted by p s = {p\,p^), can be stably supported. An input queue is said to be 
saturated if, after a head-of-the-line (HOL) packet is transmitted from the queue, there is always a packet queued 
behind it waiting to take the HOL position, i.e., the input buffer is never empty. Since each node transmits with 
probability whenever its battery is non-empty and each transmission consumes one chunk of energy, the content 
size of the battery, Bi(n), forms a decoupled discrete-time M/M/l queue with input rate 5{ and service rate pj. 
Consequently, the probability that the battery is non-empty is given by min(| i , 1). The probability of success seen 
by node i is equal to Piqi\u\ if its battery is non-empty while the battery of node i) is empty, and equal to 
Pi(l —pj)qi\{i} + PiPjqi\{i t 2} if both batteries are non-empty. Hence, the saturated throughput of node i is equal to 
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the average service rate given by 

lA = \ PiQi\{i} 



1 — min — , 1 



+ [Pi(l - Pj)Qi\{i} + PiPjQi\{i,2}} min lj | min lj 



Pj 

= min(5 i ,pj)(g i | {i} - Aixmn(6j,pj)) (13) 
for i,j G {1,2} and i ^ j. 

Lemma 4.1: The system is stable under the slotted ALOHA if the arrival rate vector A = (Ai, A2) is componen- 
twise less than or equal to the saturated throughput vector fi s = 

Proof: It suffices to show that there is no A ^ /x s , where denotes componentwise inequality, that makes 
the system unstable. The intuition is that the behavior of a node with unstable queue is statistically identical to 
that with saturated queue as time goes to infinity. This is because a queue being unstable, or equivalently transient, 
implies that its size grows to 00 without emptying with a nonzero probability. Therefore, if one queue, say queue 
i, becomes unstable at some A', then the corresponding input rate component A^ is greater than the service rate 
seen in the saturated throughput 

We formally support this argument as follows. Let be any real number satisfying < e, < /if and let 
= /■*! ~~ e «> V* £ {1,2}, such that A ^ fi s and suppose that the system is unstable with the chosen input rate 
vector. The instability of a system implies that at least one queue in the system is unstable. Let us first suppose that 
both queues are unstable such that their queue sizes, Qi(n), Vi € { 1, 2 }, grow to infinity without emptying with 
nonzero probability. Since the number of recurrences of the empty-state is finite with probability 1 , the probability 
that each battery is non-empty approaches min(^-, 1) as n — V 00, for i G { 1, 2 }. Then, the limiting expectation on 
the actual rate serviced out of queue i, denoted by limn.^oo E[fii(n)], is equal to [if in ( [TBI . On the other hand, 
from the queueing dynamics £0), for any n > 0, we have 

ra—l n— 1 

Qi{n)-Qi(0)=J2Mk)-J2^ k ) (14) 

fc=0 fc=0 

By taking expectations, diving by n, and taking a limit as n — > 00 we have 

lim E[Qi(n)} =x _ nm 1 y E[ ^ k)] (15) 

n— >oo n n— >oo n — J 

k=0 

where we use the fact that the effect of the initial queue size disappears as n — > 00 because of the stationarity of 
the arrival process. The instability of queue i implies that the left-hand-side (LHS) of (fT31 ) is strictly positive and, 
then, it must be the case that 

Ai = m! -ti> lim E[fj,i(n)] = /xf (16) 

n— >oo 

which is impossible because is non-negative. 

Consider now the case when only one of the queue in the system, say node i, is unstable and the other queue at 
node i) is stable, i.e., the limiting distribution exists only for the queue at node j while the queue at node i is 
transient. At the steady-state of queue j, the probabilities that the queue and the battery at node i are non-empty are 
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1 and min(| i , 1), respectively. These do not depend on j. In other words, node i randomly interferes with node j 
with some probability which only has the effect of lowering the average success probability of node j. Since queue 
j is stable, the average input rate is equal to the average output rate, which can be deduced by setting the LHS 
of ( fT51 ) to zero. Additional care should be taken here, because, not all the random variables with a well-defined 
distribution function have finite expectation, when the distribution is heavy-tailed with tail exponent that is less than 
1 Q3D . However, in most standard queueing systems including MIMI1, M/G/l, and G/M/l systems, the queue 
size exhibits exponential tail ||20l , ll2TTl . The queue j under consideration is a decoupled discrete-time system that 
is a variant of M/M/l system whose service is paused when the battery is empty. Since the inter-arrival times of 
energy chunks follow geometric distribution, and since the energy process is independent of the data process and 
the channel, the limiting distribution of queue j is also not heavy-tailed. Let us now compute ^j| act i ve which is 
defined as the expected number of packets that is successfully serviced from node j given that node j is active, 
i.e., given that both its data queue and its battery are non-empty. Since the probability of success seen by node j is 
PjQj\{j}, if node i's battery is empty, and pj(l — Pi)Qj\{j} + Wj<2j|{i,2}> if n °de z's battery is non-empty, we obtain 



f^j | active Pj^jlij} 



1^1 

1 — mm — , 1 

Pi 



+ [Pji 1 ~ Pihj\{j} + PiPjQj\{i,2}} mm f^-, l^j 



= Pj(lj\{j} ~ A j ram(5i,pi)) (17) 

Since the expected number of packets arriving into the data queue per slot is Xj, it follows from the property of 
a stable system, that the input rate is equal to the output rate so that the probability that node j is active at any 
given time slot is given by 

Pr[B j ^0,Q j ^0) = ^— (18) 

Afj | active 

which does not depend on 5j. As expectedly, however, if 5j is decreased, the non-active slots will be more likely 
to occur due to the emptiness of the battery, rather than that of the queue. Nevertheless, what it says is that the 
ratio between the active and non-active slots must remain the same for a given input rate Xj as long as the queue 
is stable. By noting that node j transmits with probability pj only when it is active and node i is unstable, the 
limiting expectation on the actual rate serviced out of queue i is obtained by 

lim E[fii(n)} = \ piqi\{i } ( 1 — — ) + [Pi(l - Pj)qi\{i} + PiPjqi\{i,2}] — — — j min ( — , 1 

n ^°° I \ Mj [active / M? | active J \Pi 

= mm(5 h pi) q il{i} — 

(a) 

> mm(5i,pi)(qi\{i} ~ mm(Sj,pj)) 

= Ml (19) 

where, for (a), we replaced Aj with /x* — ej and used the fact that ej is non-negative. Again, the instability of the 
queue i implies that the LHS of (fT5b is strictly positive and, then, it must be the case that 

Xi = m| - €i > lim E\pn(n)\ > n\ (20) 
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which is impossible because is non-negative. 



B. Necessary Condition 

The necessary condition for the stability of the considered system is derived through the construction of a 
hypothetical system; this hypothetical system operates as follows: i) the packet and energy chunk arrivals at each 
node occur at exactly the same instants as in the original system, ii) the coin tosses that determine transmission 
attempts at each node have exactly the same outcomes in both systems, iii) however, one of the nodes in the system 
continues to transmit dummy packets even when its data queue is empty but its battery is non-empty. The dummy 
packet transmission continues to consume one chunk of energy in the battery but does not contribute to throughput 
if the transmission is successful. Such a construction of a hypothetical system with dummy packet transmissions 
has been widely used to analyze systems of interacting queues and yields a sufficient and necessary condition 
for the stability JSJ, |[T0l , |[T3l , |[22l - |[24l . It uses the stochastic dominance technique; the queue sizes in the new 
system are path-wise never smaller than their counterparts in the original system, provided the queues start with 
identical initial conditions in both systems. However, in the case of a system with batteries, as is considered in 
this work, there exist sample-paths on which this strict path-wise dominance is violated. This is because dummy 
packet transmissions alter the dynamics of the batteries through unproductive use of their contents. For example, 
there are instants when a node is no more able to transmit in the hypothetical system due to the lack of energy 
while it is able to transmit in the original system. Being not able to transmit may imply a better chance of success 
for the other node, if the latter attempts to transmit at those instants, which causes a collapse of the sample -path 
dominance. Instead, here we use the hypothetical system of transmitting dummy packets only to derive a necessary 
condition for the stability of the original system. 

Let us define 

IZi = [a : Xi < mm(5i,pi) ( _ a ^ is W ' A -? - min (*i'Pj'X%1{i} ~ A i min (<^Pi))} ( 21 ) 

where i ^ j and i, j G { 1,2 }. 

Lemma 4.2: If the system is stable under the slotted ALOHA, then A G Ujg{i 2} 

Proof: Let us consider a hypothetical system in which node i transmits dummy packets when its packet queue 
is empty and node i) operates as in the original system, where i,j G { 1, 2 }. As mentioned earlier, all other 
random events including the packet arrivals, energy harvesting, and the decisions for transmissions have the same 
realizations as in the original system. In the hypothetical system, node i transmits with probability p, regardless of 
the emptiness of its data queue, provided its energy queue is non-empty, and each transmission consumes one chunk 
of energy. Therefore, Bi{n) forms a decoupled discrete-time M/M/l queue whose probability of non-emptiness is 
given by min(| i , 1) and node i behaves independently from node j, i.e., node i only has the effect of lowering the 
success probability of node j in the average sense. The saturated throughput of node j, therefore, can be computed 
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separately as in ( [TBI and the queue at node j is stable if 

Xj < mm(5 j ,pj)(q j \ {j ' i - Aj mm(8i,pi)) (22) 
which follows by applying Lemma |4~T1 to a single-node case. For Xj satisfying (l22l . the probability that node j is 



active is obtained as in (fT8l ) by noting that the probability that battery i is non-empty is given by min(| L , 1). Thus, 
the queue at node i is stable if 



Aj < \piqi\{i} (l — — ) + bt(! -Pj)Qi\{i} + PiPjQi\{l,2}] — — — }min(— , 1 j 

I \ P'j | active/ Mj | active J \Pi / 

/ A,A 7 - \ 
= mm(5i,pi) qi\ {i} - — r (23) 

The pair of equations (1221 and d23l ) describes the stability condition for the hypothetical system in which node 
i transmits dummy packets, which is a necessary condition for the stability of the original system for the range 
of values of Xj specified in Eq. (l22l . The reason is this: if for some Aj, queue % is unstable in the hypothetical 
system, i.e., d23l does not hold, then Qi(n) approaches infinity almost surely. Note that as long as queue i does 
not empty, the behavior of the hypothetical system and the original system are identical, provided they start from 
the same initial conditions, since dummy packets will never have to be used. A sample-path that goes to infinity 
without visiting the empty state, which is a feasible one for a queue that is unstable, will be identical for both the 
hypothetical and the original systems. Therefore, the instability of the hypothetical system implies the instability 
of the original system. ■ 

C. Proof of Theorem \3.1\ 

Here we first compute the closure of the outer bound of &(d,p) over all feasible transmission probability vectors 
pG [0, l] 2 . Therefore, any rate vector that is outside the closure is not attainable. After that, it is proven that the 
entire interior of the closure can be achieved by showing that the closure of the inner bound is identical with that 
of the outer bound. 

Note that the description on the outer bound of &(S,p) in Lemma |4~2l does not depend on S for p ^ S and 

also note that increasing pi over 5i has no effect since the value of mm(6i,pi) is bounded below by <5j. For the 

subregion IZi, i S {1,2}, let us consider the following boundary optimization problem in which we maximize the 

boundary of Aj, denoted by Aj, for a given value of Aj (j^i) as p varies[], that is 

/ AjAj \ 
max Xi = pi[ qi\ {i x J - (24) 

subject to Xj < pj(q jl{j} - AjPi) (25) 
Pi <6i,Vie {1,2} (26) 



3 Note that optimizing the boundary of a region over p is equivalent to take the closure of the region over p. 
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To maximize Aj over p, we need to understand their relationship. Note that Aj depends only on pj. Differentiating 
Aj with respect to pj gives 

dAj _ g A '9j|Q} A J (27) 

dp* * IW (^| {i } - Aj-pi) 2 

and by differentiating once again, we have 

2 \; _ 2A i A j q Mj} X j 



(28) 



Since qj\ij\ > Aj, the second derivative is negative and, thus, Aj is a concave function of pj. Equating the first 
derivative to zero gives the maximizing p* as 



and the corresponding maximum function value is obtained by substituting d29l ) into d24l ). thus yielding 



t* _/ / AjAj \ / Qi\{i}Qj\{j} ~ V A i(li\{i}qj\{j} x j 

^i, curve \ \ II a 



(30) 



Suppose now that the maximum occurs at a strictly interior point of the feasible region, i.e., p* G (0, 5i), which 
corresponds to the condition 

QmilMj} ~ A A) 2 < A c gj|{i}gj|{j} (31) 
Ai^ 10} J A, 

which is obtained by rearranging Eq. d29l ) and substituting the extreme values of p*. On the other hand, the constraint 
(|25T ) should also be satisfied for the derived p*. Hence, by substituting ( f29T > into (|25T ) and using pj < <5j, we obtain 

A . < A mM (32 ) 

?t|{i} 

Consequently, A* curve in (|30l is valid only for the range of values of Aj that satisfy both ( f3TT > and 02j ). The 
intersection of the ranges of values of Aj determined by Eqs. (I3TT ) and (l32l would be identical with the range 
specified by (fJTJ) if 5j > which is impossible because q^j > Aj while 6j < 1. Thus, if ^ > 1, where \& is 

defined in Section JIIJ the intersection is given by 



A- ■- - J 



< A < (33) 



!/} \>\ 
Otherwise, if \& < 1, the intersection is an empty set. 

Next suppose that either p* = or p* = 5i, which is the case when Aj lies outside of the range of Eq. (fJTJ). If 

Aj is on the right-hand side of the range, i.e., if Aj > <? ' I{ '^ 3|{:,> , Aj is a non-increasing function of pi since its first 

derivative in (|27T ) is non-positive. Therefore, p* = and A* = 0. On the other hand, if Aj is on the left-hand side 

of the range of (fJTJ), Aj is a non-decreasing function of pi and, hence, p* = <5j and the corresponding maximum 

function value is obtained as 

A!,nne = St (q m ^Vr) (34) 

for 

A, < q ^ q ^~ A ^ )2 (35) 
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On the other hand, the constraint d25l ) at p* = Si becomes 

Xj < 6j(qj\{j} ~ AjSi) (36) 



Thus, A* lino is valid for the range of values of Aj specified as the intersection of the ranges given by (1351 ) and (1361 1. 
If * > 1, the intersection coincides with (l33T > and, if ^ < 1, it coincides with (l36l ). To sum up, A* is obtained as 
follows: 



If * > 1, 



A* 



If * < 1, 



<\,line> Ior A J S A i9j - Kj - } 



K = X* Me , forA^^-j-A^) (38) 



Substituting i G { 1, 2 } into the above yields the description for A^ and A2 which lead us to the description for the 
stability region given in Theorem 13.11 Specifically, when \& > 1, the end points of A£ curve and Ag i ine and those of 
A* curve and A^ line meet at Pb 1 and Pg 2 (which are defined in Section HTT1 ). respectively. Furthermore, AJ curv0 and 
A2 curve are functions that are inverse of each other and they can be identically rearranged to coincide with (fT2l . 
These segments together with the axes form a closed region in the two-dimensional Euclidean space as shown in 



Fig. |2(a)| If \t < 1, the end points of AJ linc and Ag line meet at Pb 3 and, likewise, they define a closed region as 



shown in Fig. 2(b) 



What is left to be shown is the achievability of the specified region. From Lemma 14.11 we know that A < fx s 
can be stably supported. For some j e { 1, 2 }, if p ^ S, fij is written as 

= Pj(.Qj\{j} ~ A jPi) (39) 

from which we derive 

Pj = h~ ^ 

By substituting pj into the expression for /zf (i / j), we have 

" ? =»H-^A^J <41) 
which turns out to be identical to the expression for the outer boundary in (f24b by replacing /if with A, and //| 
with Aj. In other words, the operating point of the saturated system can be controlled to any point on the boundary 
of IZi by adjusting pj according to (l40b for i £ {1,2} and j 7^ i. This implies that the outer bound described by 
Lemma 14.21 can be indeed achieved, which proves the achievability of Theorem 13.11 
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Fig. 5. Two-node stability region & c (d) with finite capacity batteries where S = (0.8, 0.6) and c = (3, 3) (the case when A™ ax + A!f ax > 1) 

V. The Impact of Finite Capacity Batteries 

In this section, we consider the case where the capacity of the batteries is finite and study the impact of that on 
the previously obtained stability region. Denote by a the capacity of the battery at node i. Then, the number of 
energy chunks stored in the battery evolves according to 



i.e., the harvested energy chunks now can be stored only if the corresponding battery is not fully recharged. Since 
most of the analysis overlaps with the case of infinite capacity batteries, the result is demonstrated only for the 
collision channel model for brevity and to simplify the exposition. It becomes clear from the analysis that the 
channel with MPR capability can be handled similarly. Denote by Cj the capacity of the battery at node % and let 



Theorem 5.1: For a given energy harvesting rate 8, the two-node stability region 6 C (S) of the collision channel 
random access with batteries that have finite capacities, denoted by c = (cj, i 6 { 1, 2 }), is as follows. If A max + 
A™ ax > 1, the boundary of the stability region is composed of three segments: (i) a line segment connecting 
(0,AJj MX ) and ((1 - A!f ax ) 2 , (Af 13 *) 2 ), (ii) the curve y^i + ^% = 1 from the point with coordinates ((1 - 
A™ ax ) 2 , (A™ ax ) 2 ) to the point with coordinates ((A max ) 2 , (1 - A™ ax ) 2 ), and (hi) a line segment connecting points 
((Af ax ) 2 , (1 - Af ax ) 2 ) and (Af ax , 0). If Af ax + A!f ax < 1, the boundary is described by two straight line segments, 
namely, (i) the one connecting points (0, A™ ax ) and (A™ ax (l - A!, nax ), A!, nax (l - A™ ax )) and (ii) the one connecting 
(Ar(l - A!f ax ), A^ ax (l - A5" ax )) and (A^ ax , 0). 



Bi{n + 1) = min (Bi(n) - lj(n) + Hi(n),a) 



(42) 



max A 



5i (l-^)/(l-^). 
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Proof: We begin by noting that when a node, say node i, transmits with probability pi whenever its battery is 
non-empty in case it is either saturated or transmits dummy packets if its queue is empty, the size of the battery, 
Bi(n), follows a decoupled discrete-time M/MJl/ci model whose probability of being non-empty is given by GUI 

f (VftKi- (ft/ft)*) 



fi 



, . if Si ^ Pi 
if Si = pi 



K ci + r 

By following similar steps as in the previous section, the outer bound of 6 c (S,p) for the case of the collision 
channel model is obtained as the union of regions described by 

TU = {a : Xi < pifi (l - YZJj) ' A i ^ PiM 1 ~ (44) 
for i G {1,2} and % / j. Similarly, the saturated throughput vector of the system, which corresponds to an inner 
bound of the stability region, is obtained as /j s = (fif,i G {1,2}) with \i\ = Pifi(l — Pjfj)- By substituting 
pj = yi s -/{fj(l — pifi)) into the expression for fif, we observe that the saturated throughput vector /x s can achieve 
any point on the boundary of TZi, i.e., the rate region inside the outer bound of 6 c (S,p) can be indeed achieved. 
Let us now compute the closure of the stability region. Observe that the function p^fi is an increasing function 
of pi and, thus, its maximum occurs at p* = 1 and the corresponding maximum function value is obtained as 
Si(l — <5?*)/(l — <5j Ci+1 ), which we denoted by A™ ax . Consequently, any non-negative function value that is less 
than or equal to Af iax can be attained by appropriately selecting a value for pi between zero and one. Based on this 
observation, replace pj/j in (l44l with another variable p\ G [0, Af 1 *"]. The rest of the proof follows that of Theorem 
13.11 Specifically, we set up and solve an optimization problem similar to that of Section IIV-CI from (1241 to (|26l ) 
with <jj|{j} = 1, ^1(1,2} = 0, and pi £ [0, A? 18 *], Vi G {1,2}. Indeed, the result for the channel with the MPR 
capability can be obtained by solving the same problem as in Section [IV-CI but with pi G [0, Af ax ], Vi G { 1,2}. 
The remainder of the proof is omitted for brevity. ■ 

An example two-node stability region & c (5) with finite capacity batteries is depicted in Fig. [5] which is the 
region below the solid line. For comparison's sake, the stability region with infinite capacity batteries for the 
same parameter values is also depicted in the figure with the dotted line. The difference between the two regions, 
therefore, represents the loss due to the finite capacity of the batteries. 

Corollary 5.1: Denote by & co (5) the stability region of the slotted ALOHA with infinite capacity batteries. 
Then, for any finite capacity batteries, the relation 6 C (S) C &°°(S) holds. 

As observed in the proofs of Theorem 13. II and Theorem 15.11 we can compute (5°°(<5) and & c (5) through a closure 
operation over p varying in the rectangles [0,Si] x [0, #2] an d [0, Af 18 *] x [0, A™ 5 ™], respectively. The corollary 
follows from the fact that A™ ax is strictly less than Si for any finite 

VI. Concluding Remarks 

We studied the effect of stochastic energy harvesting, which imposes energy availability constraint on each node, 
on the stability of the slotted ALOHA. An exact characterization of the stability region was carried out for the 
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two-node case under the generalized wireless channel model with MPR capability. By comparing to the case of 
unlimited energy for transmissions, we identified the loss in terms of the size of the stability region for either 
infinite or finite battery capacities. To extend the analysis to more general networks presents serious difficulties 
of tractability due to the complex interaction between the nodes and may require approximations or alternative 
approaches that go beyond the scope of this paper. 

Appendix A 
On the Convexity of the Stability Region 

In this section, the MPR probabilities are computed for a Rayleigh fading environment and the criterion on the 
convexity of the stability region given by ^ ^ 1 is verified in terms of the values of physical layer parameters. 
The use of matched filters was implicitly assumed in using Eq. ([3]) for decoding the signal at the receiver, which 
fundamentally treats interference as white Gaussian noise. Although, techniques such as the successive interference 
cancellation |25l can improve the accuracy of the MPR capability description, comparing different physical layer 
techniques is outside the scope of our work here. 

We begin by describing the SINR of the signal transmitted from node % at the receiver as 

p 

1i\M = JT-y^ p— (45) 

where Ai is the set of nodes transmitting simultaneously, N is the background noise power, and P rx j is the received 
signal power from node i at the receiver which is modeled by 

P rX:i = ipfKr-»P tX:i (46) 

where tpi is a Rayleigh random variable with E{ipf] = 1, K is a constant, v is the propagation loss exponent, 
ri is the distance between node i and the receiver, and P tx ,i is the transmission power of node i. Let be the 
probability density function of tpf which is exponential with unit mean [17]. Then, the success probability of a 
transmission by node i when it transmits alone is computed by 



Pr 



N 



> 



(47) 

ftf(w)du (48) 

J v ± tX,2 _ 

6Nr v \ 

exp I — 77Tt~ (49) 



POO 


6Nr? 
to > — 


/ Pr 


/o 


[ ~ KP tX:i \ 



KP tx ,is 

Similarly, the success probability of a transmission by node i when it transmits along with the other node j is given 
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by 



Pr 



> 



oo roo 



Pr 



o Jo 

oo 

exp 



9(N + cujKr-vPt 



j -'txj/ 



Kr7 v Pt, 



9(N + UjKrJ v P t 



txj, 







1 + 



.Pi 



P- 



Kri u P tx ,i 

u\ -1 

exp 



9Nr\ 

kpT 



(50) 
(51) 
(52) 
(53) 



where i,j G {1,2}, j / i, and ^ and tpj were assumed mutually independent. 

From Corollary 13.11 we know that the convexity of the stability region &(S) is determined by ^ ^ 1 where ^ 
was defined as 



9l|{l} 9 2 |{2} 

By substituting the obtained MPR probabilities above into the definition of Aj, we obtain 

u\ -1 



A,; = 9- 



.P 



tx,j 



Ptx,i V "*i 



1 + 



.Pi 



tx,J 



Pt, 



exp 



Then, by substituting Aj, Vi G {1,2}, into (l54b . we express ^ in terms of the physical layer parameters as 

0Ptx, 2 r^ 2 0Ptx,ir^i 



+ 



^1 



(54) 



(55) 



(56) 



Ptx,i^2 + OPt,^ Ptx, 2 rr + 0Ptx,l*2 
Interestingly, in the case of having unlimited energy for transmissions, i.e., Si = 1, \/i G { 1, 2 }, the above criterion 
(l56l) simplifies to 9 ^ 1, which does not depend on any other parameters than the threshold 9. 
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